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We prove that all inseparable Gaussian states of two modes can be distilled into max- 
imally entangled pure states by local operations. Using this result we show that a 
bipartite Gaussian state of arbitrarily many modes can be distilled if and only if its 
partial transpose is not positive. 
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The existence of pure entangled states of two or more systems entails the possibility of find- 
ing new applications of Quantum Mechanics, in particular in the fields of computation and 
communication [1]. In practice, however, systems are exposed to interactions with the envi- 
ronment, that transform pure into mixed states, which may no longer be useful for quantum 
communication. Fortunately, there exist methods to recover pure entangled states from mixed 
ones in certain situations. These processes are called entanglement distillation (or purifica- 
tion) [2], and consist of local operations and classical communication transforming several 
copies of a mixed entangled state into (approximately) pure entangled states which can then 
be used for quantum communication. In fact, applying this method in the appropriate way 
one can construct quantum repeaters [3] that should allow efficient quantum communication 
over arbitrarily long distances even via a noisy channel. 

For this reason it is important to determine whether a given state is distillable or not. In 
general, the answer to this question is not known. At the moment we only have conditions 
that are necessary or sufficient for distillability, but not both. Clearly, only inseparable states 
can be distilled. Moreover, as shown by Horodecki et al. [4], there exists a stronger necessary 
condition, namely that p must have non-positive partial transpose (npt). In fact, there are 
entangled states which are not distillable since their density matrices remain positive under 
partial transposition [5]. Furthermore, there is evidence that this condition is not sufficient, 
since there exist npt states that nevertheless seem to be undistillable [6]. The existence 
of undistillable npt states would have interesting consequences such as non-additivity and 
non-convexity of the entanglement of formation [7] . 

On the other hand, a useful sufficient criterion, the so-called reduction criterion [8], has 
been established. It states that, given a state p on the composite Hilbert space Ti = TCa^'Hb, 
if there exists a vector G Ti. such that 

(^|trsp® 1 -pIV-) < 0. (1) 

then the state p is distillable. Here, tr^ stands for the partial trace with respect to the second 



*email: geza.giedke@uibk.ac.at 



1 



2 Distillability Criterion for Gaussian States 



subsystem. An important aspect of this criterion is that if one can find a state satisfying 
(1), then one can exphcitly construct a protocol to distill p. 

Up until now, nearly all work on the distillability problem has considered states of finite 
dimensional systems, see [9] for a current overview. In particular it was shown that states 
systems consisting of one qubit and an AT-level, N > 2 system are distillable if and only if 
(iff) they are npt [10, 6]. An alternative setting for quantum information processing, which 
considers infinite dimensional systems [continuous variables (CV) or "modes"] in Gaussian 
states is receiving increasing attention recently [11, 12]. For CV systems some distillation 
protocols for particular states have been proposed [13], and the existence of bound entangled 
states has been proved [14, 15], but the question of distillability in general has not been 
addressed. 

In this article we answer this question completely for all Gaussian states. We will prove 
that 

Theorem 1 (Distillability Criterion) 
A Gaussian state of N x M modes is distillable if and only if its partial transpose is negative. 

This shows that there are no npt bound entangled Gaussian states and, in particular, that 
for systems of 1 x TV modes all entangled states (npt is necessary for inseparability of such 
systems [15, 16, 17]) are distillable and thus useful for quantum communication. Moreover, 
our proof, which is based in part on the reduction criterion, provides an explicit protocol that 
accomplishes distillation for all those states. After introducing the necessary notation and 
properties of Gaussian states, the remainder of the paper is devoted to the proof of Theorem 
1. 

We consider bipartite systems composed of two subsystems, A and B, which consist of 
N and M "modes" [distinguishable infinite dimensional quantum systems with Hilbert space 
L^(R)], respectively. The joint system is referred to as a "iV x M system". It is convenient 
to describe the state p of such a system by its characteristic function (e.g., [18]) 

Xix) ^ tr[pD{x)]. (2) 

Here x = {qi,pi, . . . , Qn+MjPn+m) S R2Af+2M ^ ^^^^^ vector and 

D{x) = e-'^kiikXk+PkPk)^ (3) 

where Xk and Pk are operators satisfying the canonical commutation relations {h = 1). A 
characteristic function x uniquely defines a state p^. In the following we exclusively consider 
Gaussian states, i.e. states for which x is a Gaussian function of x [19] 

Xix) = e-\-^i^-id^-^ (4) 

where 7 is the correlation matrix (CM) and d S ]]^2Ar+2M ^j^g displacement. Thus, a Gaussian 
state is fully characterized by its CM 7 and displacement d. These states are of particular 
interest, since they comprise essentially all CV states that can be prepared in the lab with 
current technology. 

A matrix 7 is the CM of a physical state ifi^ (e.g. [22]) it is strictly positive, real, symmetric 
2(iV + M) X 2(Ar + M) and satisfies 

7 > J^l-^J, (5) 
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where Jn+m = ®k=i [20] with Ji = 

Simon [17] noted that for CV states partial transposition is equivalent to the orthogonal 
transformation AB{qA,PA, qB,PB) — {qa,Pa, Qb, —ps) on phase space, i.e., the momentum 
coordinates referring to B are inverted. For a Gaussian state this means that its CM is 
changed to 7 = AbI^b and the displacement to Asd. A Gaussian state with CM 7 has 
negative partial transpose (npt) iff 7 does not satisfy Ineq. (5) [17, 15], or, equivalently, iff 

7 t J^l-'J, (6) 

where J = A^JA^J is the "partially transposed" J in which the Ji's corresponding to B's 
modes are replaced by — Ji. 

The first part of the proof of the theorem is concerned with the special case of a bipartite 

two- mode Gaussian state: N = M = 1. Any such state can be transformed into what we 
called the standard form, using local unitary operations only [16, 17]. For a state in standard 
form the displacement d = and the CM 7 has the simple form 

where 




-("0 :>-(: D- 

The local unitaries needed to achieve this form are linear Bogoliubov transformations, i.e., 
generated by Hamiltonians that are at most quadratic in the operators Xi^2,Pi,2 The four 
real parameters {na,nb,kx,kp) fully characterize a 1 x 1 Gaussian state up to local linear 
Bogoliubov transformations (LLBT). They can be easily calculated from the four LLBT- 
invariant determinants det A, det B, det C, and det 7 as follows: 

na = V det A, rib = V det B, k^kp = det C, (9a) 



{naUb - kl){nanb - kp) = det J. (9b) 

Without loss of generality we choose k^ > \kp\. We call a state symmetric, if ria = rib = n, 
or, equivalently, if det A = det B. 

Now we are prepared for the proof of Theorem 1. We state the three main steps of the 
proof in three lemmas, which we prove in the remainder of this article. 
Lemma 1 (Distillability of Symmetric 1x1 States) 

A symmetric 1x1 Gaussian state with non-positive partial transpose is distillable. 
Lemma 2 (Symmetrization of 1x1 States) 

Every 1x1 Gaussian state with non-positive partial transpose can he locally transformed into 
a symmetric npt state. 

Lemma 3 (Concentrating Inseparability in two Modes) 

Every N x M Gaussian state with non-positive partial transpose can be locally transformed 
into a 1 X 1 npt state. 
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Proof of Theorem 1: The "only if" -part of the Theorem was proven proven by the 
Horodeckis in [4]. The "if" -part is clearly implied by these three Lemmas, since by Lemma 3 
the N X M case can be reduced to the 1x1 case, and that case by Lemma 2 to the symmetric 
case. ■ 
For the proof of Lemmas 1 and 2, it is useful to re-express the conditions (5,6) for 1 x 1 
states in terms of the parameters (9). We find that 7 is CM of a physical state iff 

{nanb-kl){nanb-kl) + l > nl + nl + 2kxkp, (10a) 
riaTib-kl > 1, (10b) 

and that 7 is CM of an inseparable (or, equivalently, npt) state, iff in addition it holds that 

{nam - kl){nanb - kp) + l < nl+nl- 2kxkp. (11) 

Proof of Lemma 1: For this wc use that a state is distillablc, if there exists a pure state 
lip) such that Ineq. (1) holds. This condition was proved in [8] to be sufhcient for distillability 
of finite dimensional systems. Its extension to infinite dimensions is straightforward: and 
proved in the appendix. 

Wc show now that for any symmetric npt Gaussian state p Ineq. (1) is satisfied with j^/;) 
taken as the pure two- mode squeezed state \ip) = J2n tanh" r \nn) for sufficiently large 
r > 0. Note that is a symmetric Gaussian state in standard form. We denote its CM by 
7^) and the four parameters (9) are Ua = Ub = cosh2r, kx = —kp = sinh2r. Let 7^ denote the 
CM of p. With these choices, Ineq. (1) becomes [21] 

2 [det(7trBp + TtreV')]"'^" " 4 [det(7p + 7^-)]"'^' < 0- (12) 

In the limit of large r (keeping only the leading terms in e'") this becomes after some simple 
algebra 

{n-kx){n + kp) <1. (13) 

But Ineq. (13) is implied by the inseparability criterion for symmetric states: ii Ua = rib = n 
then Ineq. (11) simplifies to 

In^ - kxkp - 1\ < n{kx - kp). (14) 

For inseparable states we observe [17] that kxkp < 0, which together with Ineq. (10b) implies 
that the LHS of Ineq. (14) is equal to — k^kp — 1 which can be transformed to {n — kx){n + 
kp) + n{kx — kp) — 1 from which Ineq. (13) follows immediately. ■ 
Since the local operation that will be shown to achieve symmctrization involves a measure- 
ment, it is more convenient to describe the state here by its Wigner function [18]. It is related 
to the characteristic function by symplectic Fourier transformation and thus is Gaussian for 
Gaussian states. The Wigner CM is related to the (characteristic) CM by = J^^~^J- 
We denote the four LLBT-invariant parameters for the Wigner CM [which are defined as 
in (9)] by {Na, Nb, Kx, Kp). We use the following easily checked facts: just as the standard 
form of 7, the standard form of can be obtained by LLBTs. A state is symmetric iff 
Na = Nb. The conditions (10,11) can be formulated equivalently in terms of the parameters 
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{Na,Ni,,Kx,Kp). While (10a, 11) are identical for the Wigner parameters, in (10b) ">" is 
changed to "<". We refer to these conditions for the Wigner parameters as (lOW, IIW) in 

the following. 

Proof of Lemma 2: If the state is not symmetric, it means that the reduced state at 
one of the two sides has larger entropy than the other. This suggests to let a pure state 
interact with the "hotter" side to cool it down. This must be done without destroying the 
ontangicmcnt of p. We proceed as follows: p is transformed to its Wigner standard form with 
parameters {Na, Nh, K^, Kp). Now assume that iVf, < Na, i.e., B is the hotter side [24]. Take 
an ancilla mode in the vacuum state and couple it to B's mode by a beam splitter [23] with 
transmittivity cos^ 0. After a measuring the ancilla's X-operator [25] results a state p with 
Wigner CM jw of the form (7) with 

A^lf c^Na + s^D, \ 

1/ I C^Na + S^NaNh J ' 



I f Nb \ ^ I f cK 



where the abbreviations c = cos^, s = sinO,!/ = s^Nb + (?, and = NaNb — K'^^p were 
used. The condition for symmetry, det j4| = detSj, requires 

tan^g= j " ' l . (15) 

Checking (IIW) for one sees that the inequality is just multiplied by (TV^tan^ 9-\-l)~^ > 0; 
therefore the transformed state is inseparable iff the original one was inseparable. It remains 
to show that there always exists a 6* to satisfy (15), i.e., that the right hand side of Eq. (15) is 
positive. The numerator is positive since we have chosen Nb < Na, the denominator is positive 
since Ni, < Na and the second part of condition (lOW) imply that {Na — D^Nb) > and the 
first part of (lOW) assures that {Na - D^Nb){Nb - DpNa) > {NaK^ + NbKpf > 0, hence 
all Gaussian states in Wigner standard form can be symmetrized this way. But since every 
Gaussian state can be brought into Wigner standard form by local unitaries, this completes 
the proof of Lemma 2. ■ 
To finish the proofs, we now turn to the general case of x M modes. Let 7 be the CM 
of a npt state. 

Proof of Lemma 3: The condition (6) is equivalent to 7 ^ i J [15]. Hence, for every npt 
state with CM 7 there exists a vector z e (0^^^+^^ such that for some e > 

z^-y - iJ)z < -e < 0. (16) 

The idea of the proof is that 7 can be locally transformed such that at both sides all but 
one mode can be discarded, and the resulting (reduced) 1x1 state is still npt. Then it is 
distillable by Lemmas 1 and 2. 

Write z in Eq. (16) as ^ = z^^^ ® 2;^^^ with real and imaginary parts Zr^\ z^^\ {x = A,B). 
We can always find a z such that ^^^^ and ^j^^ are not skew-orthogonal, i.e. {zr''^)'^ Jz^''^ ^ 
for both x = A,B [26]. 
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Now we have to find two symplectic transformations Sx,x = A,B that map the span of 
to the span of {61,62}, where 61 = (1,0, ...,0) etc. After performing the local 
transformation S = Sa ® Sb, both A and B can discard all but the first mode of their 
systems and still have an npt entangled (and thus distillable) state. That such symplectic 
transformations Sa,Sb always exist is seen as follows: let /i = 0^/2 = —Zi/{zJJzi); we 
can always extend /i, /2 into a symplectic basis [27] fk-k = l,...,2n such that J/2i+i = 
Ski, .f2k'Jhi = = f2k+iJhi+i- Then S defined by Sck = fk is symplectic [27] and maps 
spa,TLi{zr, Zi} to span{ei,62}, i.e. 



Using Eq. (17) wc sec that only the matrix elements (j)ki with k.l = 1, 2, + 1, iV + 2 
contribute to the Ihs of Eq. (18). Thus Ineq. (18) does not change if we replace 7 by the 
two-mode CM ■jred obtained from 7 by discarding all rows and columns referring to modes 
other than 1 and N + 1. This is the CM of the state in which A and B discard all but their 
first mode each. Thus Ineq. (18) shows that the state pred corresponding to jred is npt. But 
Pred is a two-mode state and thus distillable by the first part of the proof. ■ 
Note that all the operations needed to transform a general N x M npt state into a sym- 
metric 1x1 entangled state can be implemented quantum optically with current technology: 
they require nothing but squeezers, beam splitters, phase shifters [23], homodyne measure- 
ments, and the discarding of subsystems. Once a state has been transformed to symmetric 
standard form, the protocol of Ref. [8] can be used to obtain maximally entangled states in 
a finite dimensional Hilbert space. While a practical distillation protocol for such Gaussian 
states remains to be found (see however [13]), it is worth noting, that the main part of the 
universal protocol of [8], namely the filtering operation and the joint measurement, are for 
symmetric Gaussian states implemented by the procedure of Duan et aL[13]; for details see 



In conclusion, we have answered the distillability question for all Gaussian states: such 
states arc distillable if and only if they are npt. In particular, all entangled Gaussian states 
of 1 X A'^ modes are distillable, and there exist no npt bound entangled Gaussian states. 
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Appendix A 

We show that Ineq. (1) implies distillability even for dimH = 00. Let {|fc) : = 0, 1, ...} be 
an orthonormal basis of Ti, let Tin = span{|0) , |1) , |n)}, P-h„ the orthogonal projector on 
Tin , and let p be a density matrix on 7i (g) 7i. Let £ (p) = trg p® 1 — p he the map occuring on 
the lefthand side of (1). Assume that 3\ip) £H,e > such that {tp] £{p) IV') < — e < 0. Since 
Pn = PUnPPUr^ converges in the weak topology to p, there is A'' > such that £{pn) IV') < 
— e/2 for all n > N. Thus p can be projected by local operations onto a distillable state />jv 
and is therefore itself distillable. ■ 



